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Given the phenomenological success of the Nambu-Jona-Lasinio model in 
describing the meson physics in the low energy limit, it is tempting to find 
■ the fully relativistically structured nucleon solution in the same model un- 

der the similar approximation employed in the mesonic sector. To achieve 
this goal we need to solve a relativistic Faddeev equation. The factorizabil- 
ity of the two-body T-matrix reduces the three-body Faddeev equation to a 
tractable two-body Bethe-Salpeter equation. The reduced equation is then 
solved numerically. Our result indicates that the nucleon consists of three 
loosely bound constituent quarks. 
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I. INTRODUCTION 



One of the most important feature of QCD is the chiral symmetry and its dynamical 
breaking, which is expected to dictate the low energy hadronic physics. There exist work, 
such as the QCD sum rule the instanton liquid model and an explicit lattice QCD 
simulation via cooling technique directly or indirectly confirming this expectation. The 
Lagrangian introduced by Nambu and Jona-Lasinio 0] long time ago conveniently mimics 
such an essential aspect of QCD in the low energy limit. Models based on the NJL type 
of Lagrangians have been demonstrated to be very successful in describing the low energy 
mesonic physics On the other hand, due to technical reasons, these models are much 
less effective in describing low energy physics involving baryons. It is very often that extra 
assumptions beyond these models have to be used in order to make concrete predictions in 
the baryonic sector. 

While there is very little doubt that the NJL type of models could support bound bary- 
onic states, the direct approach in solving a three-body problem has only been attempted 
recently [^] with approximations apparently quite different from that of employed in the 
mesonic sector. An important point is that approximations in the baryonic sector have to 
be consistent with chiral symmetry, for example, the nucleon solution should approximately 
satisfy the Goldberger-Treiman relation |^. Otherwise the very essence of the NJL model, 
the chiral symmetry, is ruined by the ad hoc approximations. Other indirect attempts in 
finding the nucleon solution in the NJL-like models, such as the non-topological soliton 
approach P], the bosonization approach |^ and undoubtedly others can be found in the 
literature. 

In this paper we undertake the task of finding a nucleon-like solution in the NJL type of 
models. First we derive the three-body Faddeev equation in the valence constituent quark 
approximation by ignoring the three-body irreducible graphs. Due to the heavyness of the 
constituent quark, this approximation is expected to be good at low energies, as shown 
in the mesonic sector. By observing that the two-body diquark T-matrix has a separable 
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form the Faddeev equation can be reduced to an effective two-body Bethe-Salpeter equation 
with an energy dependent interaction. Then the reduced problem is solved numerically, 
without any further approximations. Although we can not explicitly show that our solution 
respects the exact chiral symmetry, in contrast with the meson solutions in the Hartree-Fock 
approximation, we believe that our work is a step forward in the right direction. So long as 
we can find weakly bound nucleon-like state of three constituent quarks, the chiral symmetry 
should be well protected, since the chiral symmetry is exact at the constituent quark level 

i- 

This paper is organized as follows. In section 2 we first introduce the model we explicit 
consider and then briefly review the two-body sector to fix parameters in the model. In 
section 3 the derivation of the three-body Faddeev equation and its reduction to the effective 
two-body equation are presented. The numerical technique involved in solving the reduced 



fully relativistic Bethe-Salpeter equation, based on the work of Rupp and Tjon |10 



IS 



recapitulated and then applied to our case in section 4. A summary and some outlook 
follows in section 5. 



II. TWO-BODY SECTOR 

The Lagrangian we consider is the two flavored Nambu-Jona-Lasinio model given by 

L = iJij^d^^jj + G,[{4j4jf + i^tj^Ta^f] - GS-f,{>^A/2W, (2.1) 

where ip is the quark fleld, (a = 1, 2, 3) and (^4 = 1, 2, ■ ■ ■ , 8) are the generators of the 
flavor SUf{2) and color SUc{S) groups respectively. Small current quark masses are ignored 
for simplicity. Since the coupling constants Gi and G2 have negative mass dimension, this 
model is not renormalizable. An appropriate ultraviolet cutoff procedure has to be specifled 
in order to make the model well deflned. In this work we insert a form factor g{k) = g{—k), 
whose functional form will be eventually taken to be a four-momentum cutoff A in Euclidean 
space for convenience, at every fermion vertex in the loop integrals. 
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The justification of using Eq.( |2.1| ) to model tlie low energy physics of the strong in- 
teraction and its phenomenological success in mesonic channels were well studied in the 
literature. A recent review can be found in reference 0. Although our primary goal is to 
find three-body baryonic solutions in this model, it is adequate to recapitulate the essential 
features of this model in the meson sector, which is used to fix all the parameters but one 
in the model. Then the two-body T-matrix in the scalar-isoscalar diquark channel, which 
consists of an essential component of the three-body Faddeev equation, will be derived. 

A. Meson Channel 

The most important feature that makes the model resembles QCD at low energy domain 
is that the NJL model and QCD share the same chiral symmetry and its dynamical breaking. 
The manifestation of this phenomenon in the NJL model is that the massless quarks acquire 
dynamical masses through the following self-consistent gap equation, when only the fermion 
bubble chain graphs are included or in the Hartree-Fock approximation, 

where m is the constituent quark mass, which is related to the fermion condensate (V'V') by 
m = -(aid + a2G2){ip^)/{N^Nf). For = 3 and Nf = 2, = 13 and aa = 8/3. 

As a consequence of the chiral symmetry breaking the pion emerges as the massless 
Goldstone boson, which manifests itself explicitly as the massless pole in the quark-antiquark 
two-body Tqg-matrix in pseudoscalar channel. If again only the fermion bubble chain graphs 
are retained or in the RPA approximation, Tgq given by Fig. [l| can be readily calculated. 
The residue of T^^- matrix at this pole, F^, has the form 

= g^gg[-ilc i-f,], (2.3) 

where (yf^^g is the pion-quark-antiquark coupling constant. The pion decay constant, is 
defined through the axial- vector current matrix element. 
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— T 

if7zPt,Sab = (0|^/'7/.75^V'kb(p)). (2.4) 

Using the chiral Ward identity, or the Goldberger-Treiman relation at the quark level, 
fngnqq = One Can easily find 

In arriving at the above result the on-shell condition p"^ = = has been used. 

There are three parameters in the model, two couplings Gi, G2 and the cutoff A. By 
equating /,r and m or (ipip) to the phenomenological values through Eq.( p.2|) and Eq. ( p.5|) 
we can fix two of them, which we pick A and G = aiGi + (22^2. This more or less fixes the 
theory in the mesonic sector. The last parameter t] = G1/G2 is left free to vary. 



B. Diquark Channel 

If we use the same fermion bubble chain approximation in the quark-quark sector, we 
can easily calculate the corresponding T-matrix. In the color 3 scalar-isoscalar channel the 
T-matrix has the structure, when ignoring the mixing with other channels (for example, the 
color 3 vector-isoscalar channel) 

T^^'^'ip) = zR{p) [Ay ® r2 ® G^,f[X^2 ®T2® ^75]'^ (2.6) 

where A^^' = (A^ — A-^)/2, G = 27072 is the charge conjugation matrix and a, b label all the 
color, fiavor and Dirac indices. The scalar function R{p) can be obtained straightforwardly 
by summing the fermion bubble chain, yielding R{p) = G'/{1 — G'J{p)) with G' = {biGi + 
^2^2)74 and 

J{p)=^^ J -0-y{k)TT[Gj,SF{k+p/2)Gj,S^{-k + p/2)], 

where Tr denotes the trace in Dirac space and Spik) and Sp{k) are the constituent quark 
propagator and its transpose (in Dirac space) respectively. Furthermore we have hi = 4, 



Whether there exists diquark bound states in this model depends on whether R{p) de- 
velops poles in the time-like region. As shown in [|1T| it is possible by varying rj to find a 
bound diquark state in this channel. It should be emphasized that the existence of such a 
diquark bound state is not a necessary condition for the existence of a three-quark bound 
state, though it might be useful to utilize the diquark concept phenomenologically to ex- 
plain certain scaling violations in lepton-nucleon experiments. In this paper the diquark 
state is merely an intermediate device in setting up the three-body Faddeev equation. The 
phenomenological relevance of the diquark will not be pursued here. 

III. THREE-BODY SECTOR 

Given the fundamental four-fermion vertex by the Lagrangian and the quark-quark two- 
body T^g-matrix, and ignoring the three-body irreducible graphs, the three-body T-matrix 
can be solved from the Faddeev equation by iterating the fundamental vertex and the two- 
body Tqq-matrix. Throwing away the three-body irreducible graphs is in some sense equiva- 
lent to ignoring the non-valence constituent quark loops in the iteration process. Due to the 
heavyness of the constituent quark mass (300 ~ 400MeV) this approximation is justified in 
the low energy region. Of course, one should realize that we do not invoke more approxima- 
tion here. Essentially the same kind of approximation was used in the mesonic and diquark 
cases. 

A. Faddeev Equation 

Since we are only interested at the moment in the three-body bound state, we only need 
to consider the homogeneous Faddeev equation. If the full three-body amplitude F-^''^ (with 
/ and d being the external fiavor and Dirac indices) is decomposed as a sum of three partial 
amplitudes F, {i = 1, 2, 3), with 

Tf/ ^ e,,,,,3r/^^^(C75)'^^"«5^^^F« (pi,p2,P3), (3.1) 
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and similarly for r2 and Fs by cyclically permuting (1,2,3), then these partial amplitudes 
satisfy the following integral equation, 

T^^\pi,P2,P3) = '^9{pi -P2)iR{pi +P2)x 

"^'^^ 9ip[ - P2)[Ci,SUp2)Ci,SM)]T^'\p[,p'„ps) (3.2) 



(2vr) 

+ / ^^?(K-P2)[C755j(K)C755F(p^)]r(^)(K,P2,P3)}. 

The notation of the above equation is depicted in Fig. ||. The factor of 2 in Eq. (|3.2| ) arises 
from the color sum eacic2^bciC2 = 25^^. Though formally similar to the non-relativistic Fad- 
deev equation, Eq.(|3.2|) is exact within the approximation mentioned above. An analogous 



equation with scalar particles was considered by Rupp and Tjon in a different context [jTO 
Since we explicitly included the color, flavor and Dirac structures in the definition of the 
two-body Tgg- matrix and three-body amplitudes F{'2 3, the recoupling-coefficient matrix has 
already been automatically taken into account in Eq.(|3.2|). 



B. Reduction to an effective Bethe-Salpeter equation 

If the two-body T^^-matrix involved has a general form, it would be a formidable task 
to find the solution for Eq. (|3.2| ). The crucial observation is that Tgg-matrix has a factorized 
form and hence we are only dealing with the so-called separable situation. The separatability 
of the two-body interaction leads to a reduction of the three-body problem to an effective 
Bethe-Salpeter equation. As a matter of fact, this reduction has already been hinted by the 
explicit form of Eq.(|3.2|). More concretely, the three-body amplitudes can be written as 



r£(pl,P2,P3) = '^dd'{Pl)9{P2 - P3)R{P2 + Ps), (3.3) 
and similarly for F*^^''^-', with \E' satisfying 

*(P3) = 4^1 ^^7(K -P2)^(P2 +P3)^7(P^ -P3)[C'755J(p^)C'755f(K)]^(K), (3-4) 

as a matrix equation in Dirac space. When deriving the above equation the quarks are 
treated as identical particles. 



Diagrammatically, Eq.( p.4|) can be represented by Fig. ^ which looks hke a boson (with 
propagator R) couphng to a third quark to form a three-body bound state. However, this 
ought to be distinguished from identifying the "boson" as the diquark bound state. The 
reduction of the three-body Eq. (p.2|) to the effective two-body Eq.( p.4| ) does not depend on 
whether the diquark channel has a pole, but rather on the separ at ability of T^g-matrix. 

Introducing the equal-mass Jacobi momentum variables q and q', 

P P 

P3 = y - g; p'i = ^~ q] (3.5) 

where the total momentum P is given by 

P = Pi +P2 +P3 =p'i +P2 +P3; (3.6) 
we find that the reduced Bethe-Salpeter equation can be written as 

^(i^, = ^ / AV(g, g'; P)R{^P + q')Km{P, q'), (3.7) 
where we have defined an energy dependent interaction 

The Dirac structure of the kernel is contained in the operator 

K = [C75(7V2 + ^)C75(7P'i + m)]. (3.9) 
Using well-known properties of the charge conjugation operator C, this simplifies to 

K = (7P2 + "^)(7p'i + (3.10) 

In view of Eqs. ( |3.5| - |3.(j| ) the various momenta present in Eq. ( |3.7| ) can be expressed in terms 
of the Jacobi variables q, q' and total momentum P. 

C. Decomposition of the Reduced Amplitudes 

To see the Dirac structure more clearly let us reduce the operator K into the Pauli form. 
Using the p -spin notation of Ref. for the upper and lower compoenents of four-spinors, 
we get for the matrix elements K{p, p') (with p, p' = ±) 
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+) = {P20 + ^){p'io + m)- CT.P2'?-P'i (3.11) 
K{+, -) = -(P20 + - {m-p'io)^-P2 

K{-, +) = {p[q + m)a.p'2 - (m - p'2q)^-Pi 
K{-, -) = (-P20 + m){-p[o + m) - a.p^a.p^ 

The simplest approximation which can be made is to neglect the lower components, i.e. the 
kernel is replaced by K{+, +). The resulting eigenvalue equation becomes in this case 

^ j d^'Viq, q'- P)R{^P + g')[(P2o + ^)(p'io + ^) - ^■P'2^-P'Mq') = Ax(g), (3.12) 

where a physical bound state solution corresponds to the eigenvalue A = 1. Assuming we 
are in the overall three-quark cm. system P = {y/s, 0), we see that there are two classes of 
solutions to Eq. (|3.12|) 

Xi = <^>i(go,k1) (3.13) 

X2 = ff-q ^2{qo, \q\) 

These classes are not coupled to each other in the integral equations. This is due to parity 
and angular momentum conservation. Xi is a s-wave solution with (/ = 0, s = 1/2, j = 1/2) 
and X2 is a p-wave (/ = 1, s = 1/2, j = 1/2). In view of the simple form of the x's the 
angular integration in Eq. (|3.12|) can explicitly be carried out. As a result, we obtain a 
two-dimensional integral equation of the form 

Xn{q) = —J dq'J q"dq'Vn{q, q'; P)R{-P + q')Xn{q') (3.14) 

ZTT'^ J-00 JO O 

where q = \q\ and 

with X = cos(^qq') and Tr2 is the trace to be taken in Pauli space. Furthermore, the operator 
On is given by 1/2 and {a. q')/ {2(f) for n=l and 2 respectively. 
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This analysis can be extended to the full equation. From Eq. ( |3.11|) we see that the Pauli 
spin dependence in \E' can be either the unit operator or a.q. In view of parity conservation 
there are also two classes of solutions, which are given by four spinors of the form 



^1 



and 



0i(go,g) ^ 

^a.q (j)2{qo,q) J 



^ a.q MQo^Q) ^ 



(3.16) 



V 



(3.17) 



The vertex functions and \l/2 are again not coupled to each other. With this form for 
\1/ a partial wave decomposed set of coupled integral equations can be derived. Inserting 
Eq. ( gl6D in Eq. (|3) we get 



with n = 1,2 and 



27r3 



E/ dq'J q"dq'VUq,q';P)R{-P + q')<Pm{<l') (3-18) 



K™(g,g';P) 



^ , 9{P'1-P'2)9{P'2-P3) 



-1 {p'2 — m'^){p'i — m?) 
The explicit expression for the matrix Knm can be determined by noting that 

Wi = [ — - — 01 -7.g — - — <p2\w, 



(3.19) 



(3.20) 



where w is a, four-spinor with every component equal 1. The operator (1 + 70) /2 is in the 
three-quark cm. system nothing else as the projection operator Vt = {M^q + '~fP)/{2Mi\f). 
With this we can now calculate Knm by projecting out the Dirac form on Q and 7.5. In so 
doing we get for the Dirac part of the kernel 



Kim = ^Tr[ii'Km(g ')] 



(3.21) 



2m 



Tr[7.g/t2(g)] 

where ki = Q, K2 = 1-q^ ■ Eq. (|3.21|) can be evaluated in a straightforward way. We find 
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Kn = -gogo ' + qo{m + ^) - q - q.q ' + (ra + (3.22) 

K,, = -(go - 2^)q- + (go ' + m - ^)(g-'.g-") 

K2i=qo -m go + 2^^ 

3 3 g^ 

K,, = q'' + [-go(go ' - m + ^) - g- + (m - 



In a similar way the coupled set of equations can be derived for \I'2. It should be noted that 
possible solutions of this type correspond to p-wave like states and as a result are expected 
not to be the ground state of the three quark system due to the centrifugal term. Since we 
are interested in this paper in the nucleon, it is natural to confine ourself to the solutions of 
the s-wave type, given by \I'i. 

IV. CALCULATIONS 



Following Ref. |T0| the resulting integral equations can be studied by performing a Wick 
rotation of the go and q^ variables to the complex plane. Assuming that the diquark system 
supports a boundstate at Mgg, we find that at the threshold point of quark-diquark scattering 
a pinching singularity can occur in the kernel of Eq. ( |3.18D at go = go = |(2m — Mgg). It can 
readily be verified that in the triquark boundstate region, corrresponding to < niq + Mqq, 
the go and q^ variables can be rotated to a path going through the point go and parallel 
to the imaginary axis without encountering any singularities in the kernel. In so doing 
we implicitly assume that eventual singularities in the form factors g{q) do not cross the 
imaginary go axis. Furthermore the arguments of the form factors are approximated by 
neglecting the go dependence. The resulting Euclidean form of the integral equation is 
regular in the boundstate region and as a result it can in principle be solved by standard 
discretization procedures. Because of the actual size of the resulting matrix equations we 
have adopted the method described in Ref. [0. The perturbation series is determined by 
iterating the equations, while the occurring two-dimensional integrals are evaluated using 
standard Gaussian quadratures. From this series the energy of the boundstate is determined 
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using the ratio method of Malfliet and Tjon [O]. It should noted that as a byproduct also 



the corresponding wavefunction can be found in this way. 

There are several parameters in the model: the cutoff A, the coupling constants Gi and 
G2- The overall mass scale can be set by the choice of the the cutoff mass A. There are 
two constraints we would like to satisfy. From the pion decay constant = 93MeV, we 
can determine according to Eq. (|2.5| ) the constituent quark mass m. Taking a value of 
A = 750MeV we find m = 375MeV. Decreasing A for instance to 739 MeV, the quark mass 
increases to m = 400 MeV. Secondly, from the self-consistent mass gap equation, the value 
of G = 13(ji + 8/3(^2 is fixed. As a consequence the only free parameter is the ratio G1/G2, 
which can be used as the parameter to determine the diquark mass. In Fig. |] the diquark 
mass dependence on this ratio is shown. 

Once the parameters of the model have been fixed we may study the three quark bound 
state. In Table | we list the diquark masses needed to get a nucleon solution at = 939 
MeV in three approximations, non-relativistic (or static limit) Ku — > 4m^; one channel 
defined by Eq.( ^.12| ) and two-channel defined by Eq. ( |3.18| ). As one can see, a stable nucleon 
solution always requires the scalar-isoscalar diquark state lie below two-quark threshold. 
The binding of the three quark system clearly depends on the choice of the diquark energy. 
In Fig. ^ are shown for two cases of A the results of the calculated mass of the three quark 
ground state as a function of the diquark mass (solid line). Also are plotted the results when 
we only keep the s-wave components of the three quark wave function (dot-dash line) and 
the static limit of Ku Am? (dash line). From this we see that at lower diquark masses 
the static limit predicts a substantially deeper binding than the full 2-channel result and 
hence it can be an unreliable approximation. 

To have a feeling on the quality of the static and one channel approximations we list 
in Table |1| the nucleon masses with diquark mass fixed at the value where the two-channel 
calculation would yield = 939 MeV. It is clear from Table ^ that the relativistic forms 
give rise to less attraction, leading to a slightly higher lying ground state, though there is no 
qualitative difference from the static approximation From these results we may conclude 
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that over a wide range of Mqq a stable nucleon solution indeed exists in the considered NJL 
model. Increasing the diquark mass leads to a weakening of the quark-quark interaction and 
as result the nucleon mass increases. 

In the range of diquark mass we considered, the existence of the nucleon solution near its 
experimental value required about 150 to 300 MeV binding in the scalar-isoscalar diquark. 
This kind of diquark clustering is also observed in a recent instanton model calculation by 
Shuryak et al [|l4l in the nucleon channel and qualitatively confirmed by the lattice simulation 



through cooling [|15[. Since the NJL type of models are practically effective theories for these 
instanton models, the similar diquark clustering in our case may not be a mere coincidence. 



V. SUMMARY AND OUTLOOK 

We have been able to demonstrate that the NJL type of models can easily accommodate 
the nucleon-like state under similar type approximations employed in the mesonic sector. 
Although we were not able to explicitly show, from the derived Faddeev equation, that the 
solution of the nucleon state satisfies the Goldberger-Treiman relation, we indeed found that 
the nucleon be a loosely bound state of the constituent quarks. In order to have the nucleon 
solution as a true bound state a bound diquark in the scalar-isoscalar channel is necessary 
in our model. 

There are clearly some interesting questions which can be addressed in such a model. 
Using the wavefunction corresponding to the three-quark boundstate, the properties of the 
various form factors for the nucleon can in principle be studied. It is also of interest to 
examine possible Delta isobar states in the same model we have considered. The mass split- 
ting between the baryon decuplet and octet constitutes a non-trivial test of the NJL type 
of models, while the mass splittings within the same baryon multiplets are less stringent 
due to the fact that the latter splittings mainly come from quark masses. Since the domi- 
nant diquark configuration in the Delta should be vector-isovector, the resulting three-body 
bound state could be a resonance rather than a bound state. The numerical method we 
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used in the nucleon case needs to be modified if tlie Delta lies in the continuum. A more 
careful examination of the compatibility of the Faddeev equation and the chiral symmetry 
could provide useful insight on how the Goldberger-Treiman relation at the nucleon level 
is realized. Finally, also pion-nucleon and nucleon-nucleon low energy scattering processes 
can in principle be studied. It is easy to anticipate that the meson-exchange potential could 
merge between nucleons, if only the valence quark lines are included at each instance in 
the Feynman graphs. Furthermore, since the nucleon in this model is a loosely bound state 
of three constituent quarks, it is very likely that there are anomalous singularities in the 
scattering processes, which could potentially modify the one-meson exchange nuclear force 
picture in the low energy limit. The role of the anomalous singularity in the form factor for 
loosely bound states in similar models was studied recently |T6[ . 
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FIGURES 



FIG. 1. Feynman graphs for Tqq in pseudoscalar channel. 
FIG. 2. Feynman graphs for Tqq in scalar-isoscalar diquark channel. 
FIG. 3. Faddeev equation for the nucleon. 

FIG. 4. Reduced effective Bethe-Salpeter equation for the nucleon. 

FIG. 5. Scalar-isoscalar diquark mass as a function of G1/G2 for two values of constituent 
quark masses. The horizontal lines indicate the quark-quark thresholds. 

FIG. 6. Nucleon mass as a function of the diquark mass, with the dash line corresponding to 
static limit Ku Am?', dot-dash line corresponding to one channel defined by Eq. (3.12), solid 
line correspoding to full two-channel defined by Eq. (3.18) solutions respectively. The dotted line 
indicates the quark-diquark scattering threshold. 
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TABLES 

TABLE L Diquark masses needed to get a nucleon solution at M^r = 939 MeV in various 
approximations . 



m (MeV) A (MeV) Mqg^ (MeV) M^^^ (MeV) M^g^ (MeV) 

375 750 579.0 570.9 576.8 

400 739 572.8 554.7 564.7 

450 728 577.0 531.5 547.8 



^static limit. 

^one channel defined by Eq. ( p.l2| ). 
*^two-cliannel defined by Eq. (|3.18D . 

TABLE n. Comparison of predictions of the nucleon mass in various approximations. The 
value of Mqn is fixed so that M^r = 939 MeV in the full two-channel calculation. 



m (MeV) A (MeV) Mgg (MeV) Mjv^ (MeV) M^v^ (MeV) 

375 750 576.8 936.4 945.3 

400 739 564.7 928.2 950.0 

450 728 547.3 892.7 957.9 



^static limit. 

^one channel defined by Eq. ( p.l2| ). 
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